IPEV: identification of prokaryotic and eukaryotic virus-derived sequences in virome using deep learning.

BACKGROUND: The virome obtained through virus-like particle enrichment contains a mixture of prokaryotic and eukaryotic virus-derived fragments. Accurate identification and classification of these elements are crucial to understanding their roles and functions in microbial communities. However, the rapid mutation rates of viral genomes pose challenges in developing high-performance tools for classification, potentially limiting downstream analyses. FINDINGS: We present IPEV, a novel method to distinguish prokaryotic and eukaryotic viruses in viromes, with a 2-dimensional convolutional neural network combining trinucleotide pair relative distance and frequency. Cross-validation assessments of IPEV demonstrate its state-of-the-art precision, significantly improving the F1-score by approximately 22% on an independent test set compared to existing methods when query viruses share less than 30% sequence similarity with known viruses. Furthermore, IPEV outperforms other methods in accuracy on marine and gut virome samples based on annotations by sequence alignments. IPEV reduces runtime by at most 1,225 times compared to existing methods under the same computing configuration. We also utilized IPEV to analyze longitudinal samples and found that the gut virome exhibits a higher degree of temporal stability than previously observed in persistent personal viromes, providing novel insights into the resilience of the gut virome in individuals. CONCLUSIONS: IPEV is a high-performance, user-friendly tool that assists biologists in identifying and classifying prokaryotic and eukaryotic viruses within viromes. The tool is available at https://github.com/basehc/IPEV.

PADLOC: a web server for the identification of antiviral defence systems in microbial genomes.

Most bacteria and archaea possess multiple antiviral defence systems that protect against infection by phages, archaeal viruses and mobile genetic elements. Our understanding of the diversity of defence systems has increased greatly in the last few years, and many more systems likely await discovery. To identify defence-related genes, we recently developed the Prokaryotic Antiviral Defence LOCator (PADLOC) bioinformatics tool. To increase the accessibility of PADLOC, we describe here the PADLOC web server (freely available at https://padloc.otago.ac.nz), allowing users to analyse whole genomes, metagenomic contigs, plasmids, phages and archaeal viruses. The web server includes a more than 5-fold increase in defence system types detected (since the first release) and expanded functionality enabling detection of CRISPR arrays and retron ncRNAs. Here, we provide user information such as input options, description of the multiple outputs, limitations and considerations for interpretation of the results, and guidance for subsequent analyses. The PADLOC web server also houses a precomputed database of the defence systems in > 230,000 RefSeq genomes. These data reveal two taxa, Campylobacterota and Spriochaetota, with unusual defence system diversity and abundance. Overall, the PADLOC web server provides a convenient and accessible resource for the detection of antiviral defence systems.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes.

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

Prokaryotic viruses impact functional microorganisms in nutrient removal and carbon cycle in wastewater treatment plants.

As one of the largest biotechnological applications, activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses, with 10-1,000-fold higher concentrations than in natural environments. However, the compositional variation and host-connections of AS viruses remain poorly explored. Here, we report a catalogue of ~50,000 prokaryotic viruses from six WWTPs, increasing the number of described viral species of AS by 23-fold, and showing the very high viral diversity which is largely unknown (98.4-99.6% of total viral contigs). Most viral genera are represented in more than one AS system with 53 identified across all. Viral infection widely spans 8 archaeal and 58 bacterial phyla, linking viruses with aerobic/anaerobic heterotrophs, and other functional microorganisms controlling nitrogen/phosphorous removal. Notably, Mycobacterium, notorious for causing AS foaming, is associated with 402 viral genera. Our findings expand the current AS virus catalogue and provide reference for the phage treatment to control undesired microorganisms in WWTPs.

Large-scale computational discovery and analysis of virus-derived microbial nanocompartments.

Encapsulins are a class of microbial protein compartments defined by the viral HK97-fold of their capsid protein, self-assembly into icosahedral shells, and dedicated cargo loading mechanism for sequestering specific enzymes. Encapsulins are often misannotated and traditional sequence-based searches yield many false positive hits in the form of phage capsids. Here, we develop an integrated search strategy to carry out a large-scale computational analysis of prokaryotic genomes with the goal of discovering an exhaustive and curated set of all HK97-fold encapsulin-like systems. We find over 6,000 encapsulin-like systems in 31 bacterial and four archaeal phyla, including two novel encapsulin families. We formulate hypotheses about their potential biological functions and biomedical relevance, which range from natural product biosynthesis and stress resistance to carbon metabolism and anaerobic hydrogen production. An evolutionary analysis of encapsulins and related HK97-type virus families shows that they share a common ancestor, and we conclude that encapsulins likely evolved from HK97-type bacteriophages.

Prokaryotic viperins produce diverse antiviral molecules.

Viperin is an interferon-induced cellular protein that is conserved in animals1. It has previously been shown to inhibit the replication of multiple viruses by producing the ribonucleotide 3'-deoxy-3',4'-didehydro (ddh)-cytidine triphosphate (ddhCTP), which acts as a chain terminator for viral RNA polymerase2. Here we show that eukaryotic viperin originated from a clade of bacterial and archaeal proteins that protect against phage infection. Prokaryotic viperins produce a set of modified ribonucleotides that include ddhCTP, ddh-guanosine triphosphate (ddhGTP) and ddh-uridine triphosphate (ddhUTP). We further show that prokaryotic viperins protect against T7 phage infection by inhibiting viral polymerase-dependent transcription, suggesting that it has an antiviral mechanism of action similar to that of animal viperin. Our results reveal a class of potential natural antiviral compounds produced by bacterial immune systems.

Significance of Viral Activity for Regulating Heterotrophic Prokaryote Community Dynamics along a Meridional Gradient of Stratification in the Northeast Atlantic Ocean.

How microbial populations interact influences the availability and flux of organic carbon in the ocean. Understanding how these interactions vary over broad spatial scales is therefore a fundamental aim of microbial oceanography. In this study, we assessed variations in the abundances, production, virus and grazing induced mortality of heterotrophic prokaryotes during summer along a meridional gradient in stratification in the North Atlantic Ocean. Heterotrophic prokaryote abundance and activity varied with phytoplankton biomass, while the relative distribution of prokaryotic subpopulations (ratio of high nucleic acid fluorescent (HNA) and low nucleic acid fluorescent (LNA) cells) was significantly correlated to phytoplankton mortality mode (i.e., viral lysis to grazing rate ratio). Virus-mediate morality was the primary loss process regulating the heterotrophic prokaryotic communities (average 55% of the total mortality), which may be attributed to the strong top-down regulation of the bacterivorous protozoans. Host availability, encounter rate, and HNA:LNA were important factors regulating viral dynamics. Conversely, the abundance and activity of bacterivorous protozoans were largely regulated by temperature and turbulence. The ratio of total microbial mediated mortality to total available prokaryote carbon reveals that over the latitudinal gradient the heterotrophic prokaryote community gradually moved from a near steady state system regulated by high turnover in subtropical region to net heterotrophic production in the temperate region.

Diverse CRISPR-Cas Complexes Require Independent Translation of Small and Large Subunits from a Single Gene.

CRISPR-Cas adaptive immune systems provide prokaryotes with defense against viruses by degradation of specific invading nucleic acids. Despite advances in the biotechnological exploitation of select systems, multiple CRISPR-Cas types remain uncharacterized. Here, we investigated the previously uncharacterized type I-D interference complex and revealed that it is a genetic and structural hybrid with similarity to both type I and type III systems. Surprisingly, formation of the functional complex required internal in-frame translation of small subunits from within the large subunit gene. We further show that internal translation to generate small subunits is widespread across diverse type I-D, I-B, and I-C systems, which account for roughly one quarter of CRISPR-Cas systems. Our work reveals the unexpected expansion of protein coding potential from within single cas genes, which has important implications for understanding CRISPR-Cas function and evolution.

VIRIDIC-A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses.

Nucleotide-based intergenomic similarities are useful to understand how viruses are related with each other and to classify them. Here we have developed VIRIDIC, which implements the traditional algorithm used by the International Committee on Taxonomy of Viruses (ICTV), Bacterial and Archaeal Viruses Subcommittee, to calculate virus intergenomic similarities. When compared with other software, VIRIDIC gave the best agreement with the traditional algorithm, which is based on the percent identity between two genomes determined by BLASTN. Furthermore, VIRIDIC proved best at estimating the relatedness between more distantly-related phages, relatedness that other tools can significantly overestimate. In addition to the intergenomic similarities, VIRIDIC also calculates three indicators of the alignment ability to capture the relatedness between viruses: the aligned fractions for each genome in a pair and the length ratio between the two genomes. The main output of VIRIDIC is a heatmap integrating the intergenomic similarity values with information regarding the genome lengths and the aligned genome fraction. Additionally, VIRIDIC can group viruses into clusters, based on user-defined intergenomic similarity thresholds. The sensitivity of VIRIDIC is given by the BLASTN. Thus, it is able to capture relationships between viruses having in common even short genomic regions, with as low as 65% similarity. Below this similarity level, protein-based analyses should be used, as they are the best suited to capture distant relationships. VIRIDIC is available at viridic.icbm.de, both as a web-service and a stand-alone tool. It allows fast analysis of large phage genome datasets, especially in the stand-alone version, which can be run on the user's own servers and can be integrated in bioinformatics pipelines. VIRIDIC was developed having viruses of Bacteria and Archaea in mind; however, it could potentially be used for eukaryotic viruses as well, as long as they are monopartite.

Mapping CRISPR spaceromes reveals vast host-specific viromes of prokaryotes.

CRISPR arrays contain spacers, some of which are homologous to genome segments of viruses and other parasitic genetic elements and are employed as portion of guide RNAs to recognize and specifically inactivate the target genomes. However, the fraction of the spacers in sequenced CRISPR arrays that reliably match protospacer sequences in genomic databases is small, leaving the question of the origin(s) open for the great majority of the spacers. Here, we extend the spacer analysis by examining the distribution of partial matches (matching k-mers) between spacers and genomes of viruses infecting the given host as well as the host genomes themselves. The results indicate that most of the spacers originate from the host-specific viromes, whereas self-targeting is strongly selected against. However, we present evidence that the vast majority of the viruses comprising the viromes currently remain unknown although they are likely to be related to identified viruses.

Regulation of Low and High Nucleic Acid Fluorescent Heterotrophic Prokaryote Subpopulations and Links to Viral-Induced Mortality Within Natural Prokaryote-Virus Communities.

Flow cytometric analysis of marine prokaryotes routinely reveals two distinct clusters of heterotrophic cells referred to as high nucleic acid fluorescent (HNA) and low nucleic acid fluorescent (LNA) populations. Evidence suggests that these may represent physiologically and ecologically distinct prokaryote populations. According to the "kill the winner" hypothesis, viral lysis reduces the efficiency of the microbial loop by decreasing the biomass and activity of the most abundant and active members of a population (i.e., competition specialist). Thus, viral-induced mortality may vary according to the physiology of HNA and LNA cells, with implications for the marine carbon cycle. Here, the abundance and production of heterotrophic prokaryotic populations were assessed in the North Atlantic during two phases of the annual plankton cycle and related to bottom-up (i.e., organic carbon variability) and top-down processes (i.e., viral abundance and lytic production). Our results demonstrate that the relative abundance of HNA and LNA heterotrophic cells and heterotrophic prokaryote production vary according to organic carbon variability in the water column, which can be strongly influenced by the physical eddy field (i.e., type of eddy: cyclonic, anticyclonic, or no eddy). In addition, the abundance and lytic production of virus subpopulations were correlated with  the cellular production and abundance of heterotrophic HNA and LNA prokaryote communities. Our data suggest group- and activity-specific linkages between hosts and viruses (i.e., HNA-V1 and LNA-V2). Specifically, V1 had a greater contribution to total viral production (i.e., 2.6-fold higher than V2 viruses), similar to their putative host. Finally, we explore potential implications of group- and activity-specific linkages between host and virus groups on the flux of carbon through the microbial food web.

Uneven host cell growth causes lysogenic virus induction in the Baltic Sea.

In the Baltic Sea redoxcline, lysogenic viruses infecting prokaryotes have rarely been detected using the commonly used inducing agent mitomycin C. However, it is well known that not all viruses are induceable by mitomycin C and growing evidence suggests that changes in trophic conditions may trigger the induction of lysogenic viruses. We hypothesized that using antibiotics to simulate a strong change in trophic conditions for antibiotica-resistant cells due to reduced competition for resources might lead to the induction of lysogenic viruses into the lytic cycle within these cells. This hypothesis was tested by incubating prokaryotes obtained throughout the Baltic Sea redoxcline in seawater with substantially reduced numbers of viruses. We used a mixture of the protein synthesis-inhibiting antibiotics streptomycin and erythromycin to induce the desired changes in trophic conditions for resistant cells and at the same time ensuring that no progeny viruses were formed in sensitive cells. No inducible lysogenic viruses could be detected in incubations amended with mitomycin C. Yet, the presence of streptomycin and erythromycin increased virus-induced mortality of prokaryotes by 56-930% compared to controls, resulting in the induction of lysogenic viruses equivalent to 2-14% of in situ prokaryotic abundance. The results indicate the existence of a previously unrecognized induction mechanism for lysogenic viruses in the Baltic Sea redoxcline, as the mode of action distinctly differs between the used antibiotics (no virus production within affected cells) and mitomycin C (lysogenic viruses are produced within affected cells). Obtaining accurate experimental data on levels of lysogeny in prokaryotic host cells remains challenging, as relying on mitomycin C alone may severely underestimate lysogeny.

Prokaryotic SPHINX replication sequences are conserved in mammalian brain and participate in neurodegeneration.

A new class of viral mammalian Slow Progressive Hidden INfections of variable (X) latency ("SPHINX") DNAs, represented by the 1.8 and 2.4 kb nuclease-protected circular elements, were discovered in highly infectious cytoplasmic particles isolated from Creutzfeldt-Jakob Disease (CJD) and scrapie samples. These DNAs contained replication initiation sequences (REPs) with approximately 70% homology to those of environmental Acinetobacter phage. Antibodies against REP peptides from the 1.8 kb DNA highlighted a 41 kDa protein (spx) on Western blots, and in situ studies previously revealed its peripheral tissue expression, for example, in pancreatic islet cells, keratinocytes, kidney tubules, and oocytes but not pancreatic exocrine cells, alveoli, and striated muscle. To determine if spx concentrated in specific neurons and synapses, and also maintained a conserved pattern of architectural organization in mammalian brains, we evaluated mouse, rat, hamster, guinea pig (GP), and human samples. Most outstanding was the cross-species concentration of spx in huge excitatory synapses of mossy fibers and small internal granule neuron synapses, the only excitatory neuron within the cerebellum. Spx also localized to excitatory glutamate type synapses in the hippocampus, and both cerebellar and hippocampal synaptic spx was demonstrable ultrastructurally. Studies of two well-characterized models of sporadic CJD (sCJD) revealed novel spx pathology. Vacuolar loss of cerebellar synaptic complexes, thinning of the internal granule cell layer, and fibrillar spx accumulations within Purkinje neurons were prominent in sCJD GP brains. In rats, comparable spx fibrillar changes appeared in hippocampal pyramidal neurons, and they preceded prion protein misfolding. Hence, spx is an integral player in progressive neurodegeneration. The evolutionary origin, spread, and neuropathology of SPHINX 1.8 REP sequences opens another unanticipated chapter for mammalian symbiotic interactions with environmental microbes.

Taxonomic assignment of uncultivated prokaryotic virus genomes is enabled by gene-sharing networks.

Microbiomes from every environment contain a myriad of uncultivated archaeal and bacterial viruses, but studying these viruses is hampered by the lack of a universal, scalable taxonomic framework. We present vConTACT v.2.0, a network-based application utilizing whole genome gene-sharing profiles for virus taxonomy that integrates distance-based hierarchical clustering and confidence scores for all taxonomic predictions. We report near-identical (96%) replication of existing genus-level viral taxonomy assignments from the International Committee on Taxonomy of Viruses for National Center for Biotechnology Information virus RefSeq. Application of vConTACT v.2.0 to 1,364 previously unclassified viruses deposited in virus RefSeq as reference genomes produced automatic, high-confidence genus assignments for 820 of the 1,364. We applied vConTACT v.2.0 to analyze 15,280 Global Ocean Virome genome fragments and were able to provide taxonomic assignments for 31% of these data, which shows that our algorithm is scalable to very large metagenomic datasets. Our taxonomy tool can be automated and applied to metagenomes from any environment for virus classification.

Avoidance of recognition sites of restriction-modification systems is a widespread but not universal anti-restriction strategy of prokaryotic viruses.

BACKGROUND: Restriction-modification (R-M) systems protect bacteria and archaea from attacks by bacteriophages and archaeal viruses. An R-M system specifically recognizes short sites in foreign DNA and cleaves it, while such sites in the host DNA are protected by methylation. Prokaryotic viruses have developed a number of strategies to overcome this host defense. The simplest anti-restriction strategy is the elimination of recognition sites in the viral genome: no sites, no DNA cleavage. Even a decrease of the number of recognition sites can help a virus to overcome this type of host defense. Recognition site avoidance has been a known anti-restriction strategy of prokaryotic viruses for decades. However, recognition site avoidance has not been systematically studied with the currently available sequence data. We analyzed the complete genomes of almost 4000 prokaryotic viruses with known host species and more than 17,000 restriction endonucleases with known specificities in terms of recognition site avoidance. RESULTS: We observed considerable limitations of recognition site avoidance as an anti-restriction strategy. Namely, the avoidance of recognition sites is specific for dsDNA and ssDNA prokaryotic viruses. Avoidance is much more pronounced in the genomes of non-temperate bacteriophages than in the genomes of temperate ones. Avoidance is not observed for the sites of Type I and Type IIG systems and is very rarely observed for the sites of Type III systems. The vast majority of avoidance cases concern recognition sites of orthodox Type II restriction-modification systems. Even under these constraints, complete or almost complete elimination of sites is observed for approximately one-tenth of viral genomes and a significant under-representation for approximately one-fourth of them. CONCLUSIONS: Avoidance of recognition sites of restriction-modification systems is a widespread but not universal anti-restriction strategy of prokaryotic viruses.

Seasonality Drives Microbial Community Structure, Shaping both Eukaryotic and Prokaryotic Host⁻Viral Relationships in an Arctic Marine Ecosystem.

The Arctic marine environment experiences dramatic seasonal changes in light and nutrient availability. To investigate the influence of seasonality on Arctic marine virus communities, five research cruises to the west and north of Svalbard were conducted across one calendar year, collecting water from the surface to 1000 m in depth. We employed metabarcoding analysis of major capsid protein g23 and mcp genes in order to investigate T4-like myoviruses and large dsDNA viruses infecting prokaryotic and eukaryotic picophytoplankton, respectively. Microbial abundances were assessed using flow cytometry. Metabarcoding results demonstrated that seasonality was the key mediator shaping virus communities, whereas depth exerted a diversifying effect within seasonal virus assemblages. Viral diversity and virus-to-prokaryote ratios (VPRs) dropped sharply at the commencement of the spring bloom but increased across the season, ultimately achieving the highest levels during the winter season. These findings suggest that viral lysis may be an important process during the polar winter, when productivity is low. Furthermore, winter viral communities consisted of Operational Taxonomic Units (OTUs) distinct from those present during the spring-summer season. Our data provided a first insight into the diversity of viruses in a hitherto undescribed marine habitat characterized by extremes in light and productivity.

High Incidence of Lysogeny in the Oxygen Minimum Zones of the Arabian Sea (Southwest Coast of India).

Though microbial processes in the oxygen minimum zones (OMZs) of the Arabian Sea (AS) are well documented, prokaryote-virus interactions are less known. The present study was carried out to determine the potential physico-chemical factors influencing viral abundances and their life strategies (lytic and lysogenic) along the vertical gradient in the OMZ of the AS (southwest coast of India). Water samples were collected during the southwest monsoon (SWM) season in two consecutive years (2015 and 2016) from different depths, namely, the surface layer, secondary chlorophyll a maxima (~30⁻40 m), oxycline (~70⁻80 m), and hypoxic/suboxic layers (~200⁻350 m). The high viral abundances observed in oxygenated surface waters (mean ± SD = 6.1 ± 3.4 × 106 viral-like particles (VLPs) mL-1), drastically decreased with depth in the oxycline region (1.2 ± 0.5 × 106 VLPs mL-1) and hypoxic/suboxic waters (0.3 ± 0.3 × 106 VLPs mL-1). Virus to prokaryote ratio fluctuated in the mixed layer (~10) and declined significantly (p < 0.001) to 1 in the hypoxic layer. Viral production (VP) and frequency of virus infected cells (FIC) were maximum in the surface and minimum in the oxycline layer, whereas the viral lysis was undetectable in the suboxic/hypoxic layer. The detection of a high percentage of lysogeny in suboxic (48%) and oxycline zones (9⁻24%), accompanied by undetectable rates of lytic viral infection support the hypothesis that lysogeny may represent the major survival strategy for viruses in unproductive or harsh nutrient/host conditions in deoxygenated waters.

Differential impact of lytic viruses on prokaryotic morphopopulations in a tropical estuarine system (Cochin estuary, India).

Our understanding on the importance of viral lysis in the functioning of tropical estuarine ecosystem is limited. This study examines viral infection of prokaryotes and subsequent lysis of cells belonging to different morphotypes across a salinity gradient in monsoon driven estuarine ecosystem (Cochin estuary, India). High standing stock of viruses and prokaryotes accompanied by lytic infection rates in the euryhaline/mesohaline region of the estuary suggests salinity to have an influential role in driving interactions between prokaryotes and viruses. High prokaryotic mortality rates, up to 42% of prokaryote population in the pre-monsoon season is further substantiated by a high virus to prokaryote ratio (VPR), suggesting that maintenance of a high number of viruses is dependent on the most active fraction of bacterioplankton. Although myoviruses were the dominant viral morphotype (mean = 43%) throughout the study period, there was significant variation among prokaryotic morphotypes susceptible to viral infection. Among them, the viral infected short rod prokaryote morphotype with lower burst estimates (mean = 18 viruses prokaryote-1) was dominant (35%) in the dry seasons whereas a substantial increase in cocci forms (30%) infected by viruses with high burst size (mean = 31 viruses prokaryote-1) was evident during the monsoon season. Such preferential infections of prokaryotic morphopopulations with respect to seasons can have a strong and variable impact on the carbon and energy flow in this tropical ecosystem.

A decade of RNA virus metagenomics is (not) enough.

It is hard to overemphasize the role that metagenomics has had on our recent understanding of RNA virus diversity. Metagenomics in the 21st century has brought with it an explosion in the number of RNA virus species, genera, and families far exceeding that following the discovery of the microscope in the 18th century for eukaryotic life or culture media in the 19th century for bacteriology or the 20th century for virology. When the definition of success in organism discovery is measured by sequence diversity and evolutionary distance, RNA viruses win. This review explores the history of RNA virus metagenomics, reasons for the successes so far in RNA virus metagenomics, and methodological concerns. In addition, the review briefly covers clinical metagenomics and environmental metagenomics and highlights some of the critical accomplishments that have defined the fast pace of RNA virus discoveries in recent years. Slightly more than a decade in, the field is exhausted from its discoveries but knows that there is yet even more out there to be found.